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A homogeneous non-relativistic Bose gas is investigated at finite temperature using renormalization 
group methods. The phase transition is shown to be second order and the effective chemical potential 
and the effective quartic coupling vanish at the critical temperature. We obtain the critical exponent 
v — 0.73 at leading order in the derivative expansion. 
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The remarkable achievement of Bose-Einstein con- 
densation (BEC) of alkali atoms in magnetic traps Q] has 
created an enormous interest in the properties of dilute 
Bose gases. A very recent review on trapped Bose gases 
can be found in gj. The homogeneous Bose gas at zero 
temperature was extensively studied in the fifties HQ- 
The properties of this system can be calculated in the 
loop expansion which is an expansion in powers of \J pa 3 , 
where p is the density and a is the S-wave scattering 
length. The one-loop correction to the ground state en- 
ergy was calculated in 1957 by Lee and Yang ||, and the 
two-loop contribution was recently obtained by Braaten 
and Nieto ||. In Ref. j(| Haugset et al. have studied 
this system at finite temperature, and included the daisy 
and superdaisy diagrams by self-consistently solving a 
gap equation for the effective chemical potential. The in- 
clusion of these diagrams are essential in order to satisfy 
the Goldstone theorem at finite temperature. Biljsma 
and Stoof have applied the renormalization group j|] 
to study the homogeneous Bose gas at finite tempera- 
ture. Their results clearly demonstrate that the phase 
transition is second order as expected, since this system 
is in the same universality class as three-dimensional xy- 
model which is known to have a second order phase tran- 
sition. Moreover, the critical temperature increases by 
approximately 10% compared to the noninteracting Bose 
gas Q . A review summarizing our current understanding 
of homogeneous Bose gases can be found in || . 

In the present work we reconsider the non- 
relativistic homogeneous spin zero Bose gas at finite tem- 
perature. We focus in particular on the critical behaviour 
and the calculation of critical exponents using RG tech- 
niques. The work of Biljsma and Stoof only included 
the effects of two and three-particle interactions by ne- 
glecting the running of higher order vertices. In our work 
we include the RG flow of higher order particle interac- 
tions and show that the critical exponents do not con- 
verge to the expected 0(2)-model results as successive 
terms are included. 

The Bose gas can be described by an effective quan- 
tum field theory pffl . As long as the momenta p of the 



atoms are small compared to their inverse size, the in- 
teractions are effectively local and we can describe them 
with a local quantum field theory. Since we assume that 
all the atoms have the same spin we can describe them 
by a complex spin zero field: 

^ = 7y| [^1 + iH • (1) 

The symmetries are Galilean invariance and a global 
0(2)-symmetry. The Euclidean Lagrangian reads 
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Here, p is the chemical potential. The ellipses indicate 
all operators that are higher order in the number of fields 
-0 and derivatives and satisfy the 0(2)-symmetry. The 
interaction g(ip^ip) 2 represents 2 — > 2 scattering and the 
coupling constant g is proportional to the S'-wave scat- 
tering length a: 
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In the following we consider the dilute gas pa 3 <C 1, which 
implies that we only need to retain the quartic interaction 
in the bare Lagrangian Eq. (Q) |7j . We also set 2m = 1 . 

In a field theoretic language, BEC is described as 
spontaneous symmetry breaking of the 0(2)-symmetry 
and the complex field if) acquires a nonzero vacuum ex- 
pectation value v. The propagator is a 2 x 2 matrix 
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^J[e p + V'(v) + V"(v)v 2 } [e p + V'{v)}, 
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and primes denote differentiation with respect to v /2. 

The grand canonical partition function has a path 
integral representation 



(6) 



In order to implement the renormalization group we mod- 
ify the action by adding a piece containing a cutoff func- 
tion, Ru(p), to the action 



where 



dr / d 



i d x{i?^W f ■ Wip + S[ip,^}\, (7) 



S[$,il> 



tl - 



dr / d d xL 



(8) 



The function R^(p) suppresses all modes in the path in- 
tegral with momenta p less than k and larger than A, 
where k and A are our infrared and ultraviolet cutoffs, 
respectively. In the present work we choose a sharp cutoff 
such that 
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k < p < A, 
otherwise. 



(9) 



The flow equation or renormalization group equation for 
the effective action T^^^] is obtained by taking a deriva- 
tive with respect to the infrared cutoff k [11-13]. This 
equation specifies how r^^u] changes as fc is lowered: 
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The sum is over the Matsubara frequencies which take on 
the values 2nnT for bosons and the integration is over d— 
dimensional momentum space. The trace is over internal 
indices. Also note that the flow equation explicitly de- 
pends on the cutoff function R k (see also fl4|). 

The effective action can be expanded in powers of 
derivatives: 

r^.kM = / dr d d x^Uf3, k [v] + Z^ k [v]eijVid T Vj 
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To leading order in the derivative expansion the effective 
action Tp^ is (up to a volume factor) simply the effective 



potential Up.k, and and are both equal to unity. 
The RG-equation for the effective potential reads 
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L0k(v) = yj^k + U'^ + U'^] [e k + U' Ptk ] (13) 
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This equation will be derived in (jq| . The first term on 
the right hand side of Eq. (|l2|) is the T = contribution 
to the RG-equation and includes the quantum fluctua- 
tions. This equation interpolates between the bare theory 
for k — A and the physical theory at T ^ for k = 0, 
since we integrate out both quantum and thermal modes 
as we lower the cutoff. This implies that the boundary 
condition for the RG-equation is the bare potential. In 
Ref. ]ll[ renormalization group ideas have been applied 
to relativistic A0 4 -theory using the real time formalism. 
Here, one can separate the propagator into a quantum 
and a thermal part, and the infrared cutoff is imposed 
only on the thermal part of the propagator. This implies 
that the theory interpolates between the physical theory 
at T = and the physical theory at T ^ 0. Hence, the 
boundary condition of the RG-equation is the physical 
effective potential at T = Jll]] . 

We have solved Eq. ( |l2|) for d = 3 numerically for 
different values of the temperature and the result is dis- 
played in Fig. [j]. The curves clearly show a second order 
phase transition, which is expected from universality ar- 
guments. 




FIG. 1. The RG-improved effective potential U/3 t k[v] for 
different values of the temperature. The phase transition is 
clearly second order. 

The effective chemical potential jip^ as we U as the 
coupling constant (A it (defined as the discrete first and 
second derivatives of the effective potential with respect 
to v 2 /2) are displayed in Fig. ||. We see that both quan- 
tities vanish at the critical point. Moreover, g^ k goes to 
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a non-zero constant at T r . The inclusion of wavefunction 
renormalization effects turns the marginal operator g^ k 
into an irrelevant operator that diverges at the critical 
temperature [[12|. 




the operator formalism and normal ordering so the zero 
temperature part of the tadpole vanishes). 

Expanding in powers of /3(1 — p,s fe ) and introducing 
the variables 
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the RG-equations can be written as 
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FIG. 2. The effective chemical potential \ip,h and the effec- 
tive quartic coupling g^ k near the critical temperature. Both 
vanish at T c . 

In order to locate the fixed points we write the ef- 
fective potential in dimcnsionless form and make a series 
expansion around the origin: 
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and = —g^ k . Truncating the series after two terms, 
we obtain the following set of equations 



: dk 
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A similar set of equations has also been obtained in 
Ref. JtJ by considering the one-loop diagrams that con- 
tribute to the running of the different vertices (They use 



Here, e = 4 — d. We have the trivial Gaussian fixed point 
(0, 0) as well as the infinite temperature Gaussian fixed 
point (—1,0). Finally, for e > there is the infrared 
Wilson-Fisher fixed point (e/5, e/ (5S d )) §. 

One can now calulate the critical exponent v which 
is related to the correlation length £ through 



i ~ \T - T c 



(21) 



Linearizing around the fixed point, we find the eigenval- 
ues (—1.278, 1.878). This implies that the critical expo- 
nent is v = 0.532 in agreement with the result of Biljsma 
and Stoof §. 

In order to check the convergence of this expansion, 
we have repeated the calculation including more terms 
in the Taylor expansion of the effective potential. The 
critical exponent v as a function of the number of terms 
N is shown in Fig. ||. The critical exponent v fluctuates 
around the value 0.73 and should be compared to experi- 
ment ( 4 He) and the e-expansion which both give a value 
of 0.67 @. Moreover, our result agrees with that of Mor- 
ris, who considered the relativistic 0(2)-model in 3c? at 
zero temperature [ [l3"| . The fact that the critical exponent 
oscillates as a function of N reflects that the polynomial 
expansion of the effective potential breaks down near the 



critical temperature |13,f3 17 



To improve our results for the critical exponents, we 
must go beyond leading order in the derivative expansion. 

We obtain the same critical behaviour as in the 
three-dimensional 0(2)-theory as expected from univer- 
sality. This can be seen directly, without computation, 
by considering the dimensionless form of Eq. (|l2] ) : 
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FIG. 3. The critical exponent v as a function of number of 
terms N in the polynomial expansion. 



The critical potential is now found by demanding 
that 111 
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and then expanding in powers of f3ujk'- 
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U^ k = -Y^k-S d ]n[0Q k ] . (24) 



Taking the limit (3 — > |12| and ignoring the piece which 
is independent of v, this leads to 
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This is exactly the same equation as obtained by Mor- 
ris for an 0(2)-symmetric scalar theory in three dimen- 
sions to leading order in the derivative expansion |l3| |. 
Therefore, the results for the critical behaviour at lead- 
ing order in the derivative expansion will be the same as 
those obtained in the three-dimensional (9(2)-model at 
zero temperature. 

In the present letter we have explicitly demon- 
strated that the phase transition of the homogeneous 
non-relativistic Bose gas is second order and that the 
critical behaviour is the same as in the three-dimensional 
scalar 0(2)-model. Both properties are expected from 
universality arguments. 
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One natural extension of the present work is to in- 
clude wave function renormalization effects by going to 
second order in the derivative expansion and investigate 
noncritical quantities such as the critical temperature 
and the superfluid fraction in the system. 
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